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Abstract. We propose a model of feedforward (open-loop) optical 
control of two-level atom in the linearized form. This model allows to 
express the general form of solution for the atomic level populations 
via the arbitrary shapes of the control signal. Then we make numer- 
ical investigations of different shapes for the optical control signal. 
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1. Introduction 

A wide spectrum of control methods can be discovered for the quantum 
systems. Among them feedforward (open-loop) approach seems to be 
the most natural, since an applied external field, can be easily designed 
as a time-dependent function. Here we will discuss the basic, but very 
important case of two-level atomic system controlled by modulated op- 
tical field. Our choice has been motivated by developed technique for 
practical design of external field in quantum optics. 



Recently other authors studied the control of two-level atoms in the frame 
of open loop-ideology when the controlling field was known a priori. It 
allowed to get the different forms of atomic energy spectra, producing 
7r- and 7r/2-pulses [lj, including the observation of the geometric phase 
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using stimulated photon echoes [2], taking special non-constant shapes 
of external field [5] etc. 

We propose a model of feedforward control for the density matrix in the 
linearized form. We use the "semiclassical approach" of the atom-field 
interaction, when a single quantum two-level atomic system (all other 
levels are neglected) is interacting with classical electromagnetic field. 
We use the standard notation following [I], but in our model the optical 
field plays the role of a control signal u(t) for open- loop (feedforward) 
control scheme [3]. A similar case for the probability amplitudes (with- 
out decay) is described in the frame of closed- loop scheme in [BJ. The 
present model has a decay component, because it involves the effect of 
elastic collisions between atoms. 



In Section 2 we present our dynamical model with atomic level popu- 
lation decay in generalized dimensionless form and then apply the lin- 
earized control procedure for different shapes of the optical control field 
u(t). This model allows to express the general form of solution for the 
atomic level populations via the arbitrary shapes of the control signal. 
Then in Section 3 we make numerical investigations of different shapes 
for the signal u. 



2. Feedforward optical control for two-level atom 

2.1. Dynamical control model for two-level atom in classical 
optical field: We consider the quantum two-level atomic system in the 
classical optical field E(t) (see Figure. 1 in Appendix). Let \a) and \b) 
represent the upper and lower level states of the atom, i.e., they are eigen- 
states of the unperturbed part of the Hamiltonian Hq with the eigenval- 
ues: Ho\a) = tuj a \a) and Ho\b) = tko^b). 

The equations of motion for the density matrix elements are given by 

%E 

Paa = -laPaa + ~^ (pabPba^ ~ plbPab e ~) ! 
%E 

Pbb = -Ibpbb - -j- (pabPb a e tuJt - p* ab Pabe~ luJt ) ; (1) 
iE 

Pab = -labPab ^Pab{Paa ~ Pbb)e W , 
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where pb a = P a ^ Pab is the matrix element of the electric dipole moment, 
7 a and 7& are the decay constants, ^ a b = (j a + %)/2 + 7^ , j ph is a decay 
rate including elastic collisions between atoms, and u = u a — u>b is the 
atomic transition frequency. 

Let's denote p a b = \p a b\^ and 

p + = p ba e*^ + p ab e-^ t+ V ; 

p_ ee i [p ba e - p ab e-^ t+ ®] . (2) 

Using d2J) we can re-write the system (JTJ) in the real form: 

\pab\E 

Paa = -lapaa H Z P- ] 

h 

\pab\E 

Pbb = -ibPbb z — p- ; 

n 

P+ = -iabp+ + up- ; (3) 

2\pab\E 

P- = -JabP- - Up + [Paa ~ Pbb) ■ 

For further calculations we put 7 a = 76 = 7. Then 

(paa + Pbb)(t) = e-^( Paa + p bb )(0) . (4) 

The first two equations of the system ([3]) can be combined together. 
We can put: 

Paa{t) ~ Pbb{t) = e _7 *x(t) ; 

p + (t) ee e-«y(t); (5) 
P-(t) = e-«z(t). 

By substitution of ([SD in ([2]) we can eliminate the decay 7-containing 
terms. Finally, rescaling the time by to: r = oot, and denoting the di- 
mensionless control signal by u(t) = 2\p ab \E(t)/huj and e = •y p h/uj, we 
get the simplified system 

x = u ■ z ; 

y = -e-y + z; (6) 
z = —e-z — y — u-x. 

Here the dot means the derivative with respect to the new dimensionless 
time r. We remind that x G [—1,1], since (p aa — p b b) € [—1,1], and 
{paa ~ Pbb) -> as t -> 00. 
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2.2. Linearization of control: Let's suppose that we apply the lin- 
earized form of control: 

x(t) = X q {t)+u-X 1 (t) ; 

y(r) = Y^ + u-Y^t) ; (7) 
z{t) = Z (r) + u ■ Z 1 (r) . 

We will skip all the terms of the order u 2 and elder. Then substituting 
dZ|) in (JBJ), we split our system into two parts: the free (non-controlled) 
system: 

X = 0; 

Y = - e -Y + Z ; (8) 
Z = —e-Z — Y 

and the controlled part: 

u • Xi + u ■ X\ = u • Zq ; 
u ■ Y 1 + u ■ Y x = u ■ Z l ■ (9) 
u ■ Zi + u ■ Zi = — u ■ Yi — u ■ X . 

In ([9]) we omitted the decay e-terms, because the decay is supposed to 
be a slow process to compare with the control, i.e. e and u are the small 
parameters of the same order, and the linearization deals only with their 
first orders. Then from the first equation of system ([9]) we get: 



w(t)Xi(t) = / u(t')Z (t')dt' 
Jo 

and from the first equation of system ([7]), we have 

x(t) = X (t) + / u(t')Z (t')dt' . (10) 







Now we apply the initial conditions X (0), Y o (0), Z (0) to solve the 
system ([8]): 

X (t) = X (0) = x(0) ; 

Y {t) = e- eT [Y (0) cost + Z (0) sin r] ; (11) 
Z q (t) = e- eT [Z (0) cost -y (0) sin r] . 

If we denote the phase of p a b by 0', then p+(0) = 2\p a b\ cos(0' — 4>) 
and P-(0) = 2\p a b\ sin(<// — <p). We can put for the initial condition: 
4>' = (ft, then p+(0) = 2|p a6 (0)| ee 5 and P-(0) = 0. Let's demand 
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Xi(0) = Fi(0) = Zi(0) = 0. Thus, F o (0) = 5 and Z (0) = are our 
initial conditions. 

2.3. Control signal correction: If Xq(0) = — 1 (that corresponds to 
the ground level of the atom as the initial condition), then from — 1 < 
x(t)<1 and ([TDD Jt follows: 



< f u(t')Z Q (t')dt' < 2. 
Jo 



(12) 



In other words this integral should be positive and bounded. We define 
first the arbitrary non-corrected control Uq(t) and then put 



u(t) = \u (t) \ ■ signZo(r) . 



(13) 



Then the left inequality ([T2l will be satisfied automatically. The right 
part of f|T2|) can be represented by Cauchy - Schwartz inequality: 



u(t')Z {t')dt' 



< 



u\t')dt'- / Zl{t")dt" 
n Jo 



and then we demand: 

r u 2 (t')dt' ■ f zi(t")df <4 . 

Jo Jo 
Let's check the inequality ffl"4l) : 



(14) 



Zl(t")dt" 



IW)| S 



-2et" 



e sm t dt 



<S 2 

Thus, from (HI and ([15]) 

M 2 (t')rft' < 







< 



5 2 a 



-2er\ 



2e 



5 2 (1 



-2<ET^ 



(15) 



(16) 



To satisfy ([TB"]) we also have to correct the signal m. Let's suppose that 
there are two functions: an initial arbitrary uq(t) and its corrected vari- 
ant u(t) that is bounded above by the condition f[T6l) . Of course, physi- 
cally the external optical field should follow the signal u, and the initial 
uq is only a basic model to construct the behavior of the open-loop con- 
trol field. 



(> 
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Now let's define 



A(r)= J^W-^Jt— (17) 



and 

' u(t) , A(r) < ; 



B(t) , A(r) > , 
where a positive function B(t) is defined from the next equation: 



B2{t ' )dt ' = 5H1 -e-^ ' (19) 



or 

B\r 

Thus 



d 



dr 5 2 (l-e- 2 ' 



2 „-2er 



lQe z e 



6 2 n _ e -2er\2 



For small time intervals r << 1/e we have: 5 2 t in RHS f|T5|) . and 
B(r) ~2/(S-t). 

Finally by the corrections f|T3|) and f|T8|) we have: 

x(t) = -1 + 5 ■ I d£ e~ e *'| sint'l ■ \u(t')\ . (21) 
Jo 

Eq. (!2T]) solves the problem of open-loop control in linearized form. Now 
defining the control signal u(t) we restore by (I2TT) the shape of the dif- 
ference Paaif) — Pbbit)- Their sum (jlj) is known, thus, we can find p aa (t) 
and pbbit) separately. 

3. Numerical simulation of different shapes for the 

control signal 



Now we can apply the general solution of Eq.( j21j) to study the influence 
of control optical field u on the behavior of the system (jH]). 

In the case of an ideal open-loop control the behavior of x(t) is the satu- 
ration of the population at the ground level, in other words, x(t) — > 1 as 
t — > oo. Sure, not every control will satisfy this condition. 

On. Figs. 2-6 (see Appendix) we plot the different shapes of the initial u 
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(Figs. A) and corrected u (Figs.B) control signals: constant, ramp, step, 
sine wave, and repeating sequence stair. To compare their efficiency we 
check also the corresponding time derivatives dx/dt (Figs.D). 

We can see from the plots that the ramp control in our case is defi- 
nitely more effective. The speed of the saturation for x(t) is faster for 
the signals on Figs. 3, 6. 

4. Conclusion 

Finally we can conclude that our model for the open-loop control has 
several important features: 

1. It can be easily extended for the case of multi- level atomic systems by 
adding the correspondent components in the density matrix; 

2. For the two-level system it can be re- formulated in general form if 
we propose the linear approximation of control; 

3. It can be an origin of studying the behavior of controlled non-linear 
systems in quantum optics. 
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6. Appendix: Figures 
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Figure 1 . Interaction of a single two- level atom with an 
optical field. 




FIGURE 2. Constant control signal: (A) The initial signal 
u (t); (B) The corrected signal u(t); (C) x(t); (D) The 
derivative dx/dt. 
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Figure 3. Ramp control signal: (A) The initial signal 
Ua(t); (B) The corrected signal u(t); (C) x(t); (D) The 
derivative dx/dt. 
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Figure 4. Step control signal: (A) The initial signal Uo(t); 
(B) The corrected signal u(t); (C) x(t); (D) The derivative 
dx/dt. 
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Figure 5. Sine wave control signal: (A) The initial signal 
Uo(t); (B) The corrected signal u(t); (C) x(t); (D) The 
derivative dx/dt. 




Figure 6. Repeating sequence stair control signal: (A) 
The initial signal u (t); (B) The corrected signal u(t); (C) 
x(t); (D) The derivative dx/dt. 



